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Commissioner for Patents 

P.O. Box 1450 

Alexandria, VA 22313-1450 



Dear Sir/Madam: 



This Declaration is submitted to establish invention of the subject matter of the 
present patent application prior to the publication date of February 12, 2004 of U.S. Patent 
Application No. US 2004/0027363 (hereinafter referred to as "Allen"). 

1 . The person making this Declaration is inventor David C. Collins. 

2. Accompanying this Declaration is Exhibit A to establish reduction to practice of 
the subject matter of the present patent application in the United States prior to 
February 12, 2004. 

3. Exhibit A (14 pages) includes a Hewlett-Packard Company (HP) Invention 
Disclosure and accompanying attachment (hereinafter referred to collectively as 
the "Disclosure"). 

4. The Disclosure was prepared by the inventor in the United States prior to 
February 12, 2004. 

5. The Disclosure was witnessed in the United States prior to February 12, 2004. 

6. The Disclosure was submitted by the inventor to the HP Legal Department in the 
United States prior to February 12, 2004. 

7. The Disclosure was assigned HP Patent Disclosure No. 200400519. 



8. The Disclosure describes the subject matter of the present patent application — i.e., 
the subject matter of independent claims 1,12, 20, and 25 - as follows. 



Claims 


Disclosure 


Claim 1: 


The Disclosure (pp. 1-14) describes the 
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"A method of displaying an image with a 


generation of sub-frames for display at 


display device, the method comprising:" 


spatially offset positions to form a displayed 


Claim 12: 


image. 


"A system for displaying an image, the 


Wobulation, as referenced on pp. 2, 4, and 7 


system comprising:" 


of the Disclosure, refers to the display of sub- 


Claim 20: 


an image using a display device. 


"A system for generating sub-frames for 




display at spatially offset positions to 




generate the appearance of an image, the 




system comprising:" 




Claim 25: 




"A comtyirter-readable medium ^torincr 




comnuter- executable instructions which 




when executed by a computer processing 




system, cause the system to perform a 




method of generating a sub-frame image 




which comprises sub-frames for display at 




spatially offset positions to generate the 




appearance of a displayable image, 




comprising:" 




Claim 1: 


Image data is shown on p. 5 of the Disclosure 


"receiving image data for the image;" 


(a region of interest "ROI" for a given sub- 


Ulaim iz. 


frame pixel value is shown relative to the 
image data). 


"a buffer adapted to receive image data for 
the image;" 


Image data is operated on by the set of 
computer-executable instructions listed on 


Claim 20: 


pp. 9-14 of the Disclosure. 


"means for receiving image data 




corresponding to the image;" 
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Claim 25: 

"receiving image data corresponding to the 
image;" 

Claim 1 : 

"generating first and second sub-frames, 
wherein the first and the second sub-frames 
comprise a plurality of sub-frame pixel 
values and a plurality of error values, and 
wherein at least a first one of the plurality of 
sub-frame pixel values is calculated using the 
image data, at least a second one of the 
plurality of sub-frame pixel values, and at 
least one of the plurality of error values;" 

Claim 12: 

"an image processing unit configured to 
generate first and second sub-frames 
comprising a plurality of rows of sub-frame 
pixel values, wherein each of the sub-frame 
pixel values in each of the plurality of rows is 
calculated using the image data, at least one 
sub-frame pixel value from a previous one of 
the plurality of rows, and at least one error 
value;" 

Claim 20: 

"means for generating a plurality of 
rows of initial values using the image data; 

means for accessing a row of history 
values and error values generated using the 
image data; and 

means for generating a sub-frame 



On p. 4 of the Disclosure, "[b]oth of the sub- 
frames are processed together at the same 
time, and the sub-frames are intertwined." 

The diagram on p. 4 of the Disclosure 
illustrates how the sub-frames are intertwined 
in one embodiment. 

Pp. 7-8 of the Disclosure describe the 
generation of final sub-frame pixel values for 
the sub-frames using history values (e.g., 
"Previous row of subframe values" in the 
diagram on p. 8) and error values (e.g., 
"Previous subframe errors" in the diagram on 
p. 8). 

The set of computer-executable instructions 
listed on pp. 9-14 of the Disclosure performs 
a method of generating a final sub-frame 
pixel value (rv) using other final sub-frame 
pixel values (e.g., finalO_5M, final0_4M, 
fmalO_3M, fmalO_2M, fmal0_lM, and 
final0_0M,) and initial values (e.g., 
image 1_1M to imagel_4M, image2_lM to 
image2_4M, image3_lM to image3_4M, 
image4_lM to image4_4M). 
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pixel value using the row of history values 




and error values and the plurality of rows of 




initial values." 




Claim 25: 
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corresponds to a second sub-frame using the 




image data; and 




generating a second sub-frame pixel 




value using the image data, the second 




plurality of initial values, the first history 




value, and the first error value." 




Claim 1: 


Wobulation, as referenced on pp. 2, 4, and 7 


"alternating between displaying the first sub- 


of the Disclosure, refers to the display of sub- 


frame, including displaying the first one of 


frames at spatially offset positions using a 


the plurality of sub-frame pixel values, in a 


display device. 


first position and displaying the second sub- 




frame, including displaying the second one of 
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the plurality of sub-frame pixel values, in a 




second position spatially offset from the first 




position" 




Claim 12: 




"a display device adapted to alternately 




display the first sub-frame in a first position 




and the second sub-frame in a second 




position spatially offset from the first 




position." 





9. The subject matter of the present patent application was tested as follows. 

a. The set of computer-executable instructions is listed on pp. 9-14 of the 
Disclosure embodies the method of generating a final sub-frame pixel 
value described in the Disclosure. 

b. The set of computer-executable instructions listed on pp. 9-14 of the 
Disclosure was written prior to February 12, 2004. 

c. The set of computer-executable instructions listed on pp. 9-14 of the 
Disclosure was executed with at least one test image on at least one 
computer system to generate final sub-frame pixel values for sub-frames 
for the test image prior to February 12, 2004. 

d. The sub-frames generated for the test image were displayed on at least one 
display device to confirm the successful generation of sub-frames and the 
successful reproduction of the test image prior to February 12, 2004. 

10. The execution of the set of computer-executable instructions listed on pp. 9-14 of 
the Disclosure on the computer system to generate the final sub-frame pixel values 
for the sub-frames combined with the display of the sub-frames on the display 
device demonstrates actual working conditions or a realistic simulation of working 
conditions for the subject matter of the present patent application. 
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1 1 . The successful generation o^^^ successful reproduetion of the test 

image demonstrates utility beyond a probability of failure for the subject matter of 
the. present patent application, 

12. The test results for the subject matter of the p^ 

successful generation of sub-frames and successful reproduction of the test image, 
are reproducible by executing the set of computer-executable instructions listed on 
pp, 9-14 of the Disclosure on a suitable computer system with another suitable 
image to generate final sub-frame values for sub-frames forthe^image and 
displaying the sub-frames on a suitable display device to reproduce the image. 

13, As demonstrated by this Declaration* Exhibit A> and Exhibit B, . the subject matter 
of the present patent application was reduced to practice in the United States prior 
to the publication date of February 12, 2004 of Allen. 

As a person signing below, I hereby declare; that all statements made herein of my 
own fcn^ tnade on infemrafe believed 

to be (rue; and flirther that these statements were made wiih that wilMl false 

statements and the like so made ^evpunishab^^ line or imprlspnmerit, or both, under 
Section 1001 of Title 18 of the United States Code, and that such willifiil .^s^taitei^Qts .fti^y 
jeopardize the validity of the application or any patent issued thereon. 
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200400670: A Practical Implements^... 



nvent 



Innovation Number 
200400670 



D.sC^Sun^ry 




Title A Practical Implementation of the 2 Position Multipass Center Adaptive Algorithm 

Abstract The 2 position center adaptive algorithm cannot be similified in the same fashion as the 4 

position center adaptive algorithm. This adds a complication to the actual implementation of the 
algorithm in the ASIC. This disclosure shows how to overcome this complication. This is done 
by storing both the final subframe values and the error values from the previous row. 



Attachments ^ Adaptive_KerneLWIth_History_2Pos.doc [285696 bytes] 




Inventors David C Collins 



nvention Disclosures 




'■■WW'S- 



Disclosure No. 200400670 




The information contained in this document is HP CONFIDENTIAL and may not be disclosed to others without prior authorization. 
Submit this disclosure to the HP Legal Department as soon as possible. No patent protection is possible until a patent application 
is authorized, prepared, and submitted to the Government. 


j^, General Information 


Title 


A Practical Implementation of the 2 Position Multipass Center 
Adaptive Algorithm 


Abstract 


The 2 position center adaptive algorithm cannot be similified in the 
same fashion as the 4 position center adaptive algorithm. This adds a 
complication to the actual implementation of the algorithm in the 
ASIC. This disclosure shows how to overcome this complication. This 
is done by storing both the final subframe values and the error values 
from the previous row. 






Projects 




Products 




j^J Attachments 




Attachments 


\ f Adaptive__KerneLWith_History_2Pos.doc [285696 bytes] - 
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pfjy Description of Invention 



rrODiems ooivea 


i nis invenTion snows now to implement a mutiipass center adaptive 
algorithm for 2 position wobulation using the same memory 
requirements as the simplified 4 position center adaptive algorithm 


Prior Solutions 


Most of the prior solutions for 2 position wobulation have been single 
pass. My recent disclosure M A Practical Implementation of Multipass 
Adaptive Wobultion" can only be extended to the 2 position algorithm 
if the standard adaptive method is used, but the standard adaptive 
algorithm does not perform well for some types of input. 


Description 


The complete description is contained in the attatched document. For 
the implementation of the multipass adaptive algorithm, one region of 
memory is required to store the previously processed subframe row. 
The novelty of this invention is storing the error values from the 
previous processed row as well. In particular the error value that is 
stored is as follows: error = errorjeft + 2*error + error_right By storing 
tho error values in this fashion the error values and the final 
subframe values can both be stored in an interleaved fashion in the 
same region of memory. This enables the 2 position algorithm and the 
4 position algorithm to be implemented using the same amount of 
memory. The other novelty of this invention is transforming the error 
value from a signed number containing many bits to an 8 bit number. 
A simple conversion routine is defined in the attatched document; and 
another embodiment would be to use a look up table to do the 

UUI IVfcM olUf I. 


Advantages 


The primary adavantage of this algorithm is that it enables to 2 
position center adaptive wobulation algorithm to be performed using 
the same memory requirements as the 4 position algorithm. This 
invention enables one ASIC to be development that contains both 
algorithms with out any unnecessary memory. Without this invention 
the 3 pass 2 position center adaptive algorithm would require and 
additional 20% of on chip memory. 
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Adaptive Kernel With History - 2 Position 

Background 

The following explanation is assuming 2 position wobulation. Thus, two low resolution 
subframes must be generated for each frame — one for each wobulation position. Both of 
the subframes are processed together at the same time, and the subframes are intertwined. 
Thus, every second pixel will correspond to a different subframe. The following 
diagram illustrates the idea. The grey pixels are not used for two position wobulation. 



Top Left 
Subframe 



Bottom Right 
Subframe 




To compute the optimal solution, and iterative algorithm can be used. The iterative 
approach has been called the adaptive algorithm. The standard 4 position adaptive 
algorithm is described below. 
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Iterative Algorithm 
Simulation = K s * Subframes 
Error = Image - Simulation 



Origina 



image 

muianon 
Error aV g = K e * Error 
Subframes = Subframes + ax Error avg 



In a previous disclosure, I showed how the 4 position standard adaptive algorithm could 
be implemented in one pass using a small region of interest. One of the key discoveries 
was that the pixels above the pixel of interest were only needed to generate a simulated 
image (all the subframes merge together). The following diagram illustrates the region of 
interest for one pass of the 4 position standard adaptive algorithm. 
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ROT for Error Averaging 


























ROT for Tnitial Values 









ROT for Simulation 



Since the pixels above the pixel of interest are only required for the simulation ROI, the 
calculated values for the previous row are exactly the values that are needed. This same 
approach works for the 2 position standard adaptive method. The only change is that the 
simulation kernel has a value of Vi instead of l A. This change occurs because only half of 
the high resolution pixels are non-zero. Everything else for the adaptive kernel with 
history works exactly as was disclosed before. 
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Image 



Iterative Algorithm 
Simulation = K s * Subframes 
Error = Image - Simulation 
Error avg = K e * Error 
Subframes = Subframes + ax Error aV g 



The next algorithm that I considered for 4 position wobulation was the center adaptive 
algorithm. This algorithm performed better than the standard adaptive algorithm on 
single pixel lines. The center adaptive algorithm has several drawbacks in terms of 
implementation. First it includes many more computations than the standard adaptive 
algorithm. The second disadvantage is that both the error ROI and the simulation ROI for 
a given pixel extend above and below the pixel of interest. This implies that for the 
adaptive kernel with history. Both the final subframe values and error values from the 
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previous row are required for the algorithm. Thus more on chip memory would be 
required for an ASIC implementation, 

at 
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Iterative Algorithm 
Simulation = K s * Subframes 
Error = Image - Simulation 
Error avg = K e * Error 
Subframes = Subframes + a x Error aV g 



Origina 



To overcome the implementation difficulties with the center adaptive algorithm. I 
presented a simplified center adaptive algorithm. This simplified center algorithm 
performed comparably to the center adaptive algorithm, and it was less computationally 
expensive than even that standard adaptive algorithm. The simplified algorithm is given 
below. One of the key features to note is that the error averaging step has been 
eliminated. This is what allowed the adaptive kernel with history approach to work 
without the burden of additional memory. 




Origina 



Iterative Algorithm 
Simulation = K s * Subframes 
Error = Image - Simulation 
Subframes = Subframes + ax Error 
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Unfortunately these simplifications do not give satisfying results for two position 
wobulation. This approach fails because only half of the subframe values are non-zero. 
Two of the subframes don't exist at all, and this is the same as setting all their values to 
zero. 



2 Position Center Adaptive Kernel with History 

Consider again the complete center adaptive algorithm. For 2 position wobulation, 
conceptually, half of the high resolution values are zero (two of the four subframes don't 
exist at all). Thus, the simulation kernel can be reduced to the following two kernels. 
One kernel is used when the center hi-res pixel is non-zero, and the other kernel is used 
when the center pixel is zero. In addition, the error only needs to be averaged for the 
nonzero pixel location (the error only has to be feedback for subframes that exist). 

Simulation Kernel 
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Error Kernel 
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Pixejl *2 1 
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Iterative Algorithm 
Simulation = K s * Subframes 
Error = Image - Simulation 
Error avg = K e * Error 
Subframes = Subframes + a x Error avg 



Origina 



Thus, the computational complexity is only half of the complete 4 position center 
adaptive algorithm. The obstacle of the error from the previous row still remains. The 
solution is to store both the final subframe values and the error values in one row of 
image memory. For a given high resolution row, only half of the locations are used to 
store subframe values. The other half of the values are unused or set to zero. This leaves 
half of the values for storing the errors. 

The second observation is that for a given row in the hi-res image, only half of the values 
require the error to be averaged. The following diagram is an attempt to illustrate this 
point. 
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The pink error values actual contains a summation of the error value of the left and right 
pixels in addition to its own error value. In particular, each error value adheres to the 
following formula: 

Error = l*error i e ft_pixei + 2*error + l*error r i g ht_pixei 

This matches the first row of the error averaging kernel. Thus, in one row of memory we 
can store both the required error values and the subframe values. The last little detail is 
that the error value is a signed value, and it contains more bits than a single pixel. To 
accommodate this we can use a lookup table or a simple mapping. Psuedo code for one 
such simple mapping is as follows: 



temp = error_lef t+2*error+error right; 


// 


lx 2x Ix 


temp — temp/ 4; 


// 


divide by 4 


if ( temp < -127 ) temp = -127; 


// 


clip value 


if ( temp > 127 ) temp = 127; 


// 


clip value 


temp += 127; 


// 


shift to make non-zero 



Thus, the key features to implement the 2 position center adaptive algorithm using the 
adaptive kernel with history is that both the final subframe values and the previous error 
values can be stored in one row of image memory. Thus, a 3 pass 2 position center 
adaptive algorithm can still be implemented with 1 row of history and 4 rows of image 
data. 
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Appendix 

For completeness I have include some C++ code for calculating the center adaptive kernel with history. One 
thing to note is that I need to keep track of whether I am on an odd or an even pixel because I have to keep track 
of which subframe values are unused. 



unsigned char Adapt iveCenterKernel_2Pos :: Calculate 
( unsigned char final 0, 

unsigned char image 1, 

unsigned char image2, 

unsigned char image3, 

unsigned char image4, 

unsigned char isOddPixel 

) 
{ 



int temp; 












// 


isOddPixel 


= true 








// 


0 


fO 1 


0 


fO 3 


0 


fO 5 


// 


gl 0 


0 


gl 2 


0 


gi 4 


0 


// 


0 


g2 1 


6 


g2 3 


0 


q2 5 


// 


g3 0 


0 


g3 2 


0 


g3 4 


0 


// 


0 


g4 1 


0 


g4 3 


0 


g4_5 


// 


isOddPixel 


— false 






// 


fO 0 


0 


fO 2 


0 


fO 4 


0 


// 


0 


gl 1 


0 




0 


gl_5 


// 


g2_0 


0 


g2 2 


0 


q2 4 


n 


// 


u 


g3 1 


0 


g3 3 


0 


g3 5 


// 


g4 0 


0 


g4 2 


0 


g4 4 


0 



shiftValues () ; 

final 0_5M = finalO; 
image 1_5M = image!; 
image 2_5M = image2; 
image3_5M = image3; 
image 4_5M = image 4; 

.// calculate guess for column 4 

// using pixel selection - no calculations here 
guess 1_4M = image 1_4M; 
guess 2_4M = image2_4M; 
guess3_4M - image3_4M; 
guess 4_4M = image 4_4M; 

// compute sim column 4 
int siml_4 = 0; 
int sim2_4 = 0; 
int sim3_4 = 0; 
int sim4_4 = 0; 

if ( isOddPixel ) 
{ 

siml_4 = final0_3M+guess2_3M+£inal-0_5M+image2_5M; 
siml_4 += guessl_4M«2; // multiply by four 

sim2_4 = (guessl_4M«l) + (guess2_3M«l) + (image2_5M«l) + (guess3_4M«l ) ; 

sim3_4 = guess2_3M+guess4_3M+image2_5M+image4_5M; 
sim3_4 += guess3_4M«2; // multiply by four 
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} 

else 
{ 



} 



//sim4__4 » (guass3_4M«l) + (guess4___3M«l ) + <guess4__5M«I) + (guess5__4M«l) / invalid 
sim4_4 = (gues=3__4M«l) + (guess4_3M«l) + { image 4_5M«1> + (guess4__4M«l) ; 



siml_4 = (finalO_4M«l) + (guessl_3M«l ) + <imagel_5M«l ) + (guess2_4M«l ) ; 

sim2_4 = guessl_3M+guess3_3M+imagel_5M+image3_5M; 
sim2_4 += guess2_4M«2; // multiply by four 

sim3_4 = (guess2_4M«l) + (guess3_3M«l) + ( image 3_5M«1) + (guess4_4M«l) ; 

// sim4_4 = guess3__3M+guess5_3M4guess3_5M+guess5_5M; invalid 
sim4_4 = guess3_3M+guess4_3M+image3_5M+image4_5M; 
sim4_4 += guess4_4M«2; // multiply by four 



int errl_4 = (imagel_4M«3) - siml_4; // column 4 

irit err2_4 = (image2_4M«3) - sim2_4 

int err3_4 = (image3_4M«3) - sim3_4 

int err4_4 = (image4_4M«3) - sim4_4 

// compute sirn column 3 
int siml_3 = 0; // column 3 
int sim2_3 = 0; 
int sim3_3 = 0; 
int sim4_3 = 0; 

if ( !isOddPixel ) 
{ 

siml_3 = f inal0_2M+guess2_2M+final0_4M+guess2_4M; 
siml_3 +- guessl_3M«2; // multiply by four 

sim2_3 = (guessl_3M«l) + (guess2_2M«l) + (guess2_4M«l) + (guess3_3M«l) ; 

sim3_3 = guess 2_2M+guess4_2M+guess2_4M+guess4_4M; 
sim3_3 += 'guess3_3M«2 ; // multiply by four 



} 

else 
{ 



} 



sim4_3 = (guess3_3M«l) + (guess4_2M«l) + (guess4_4M«l) + (guess4_3M«l ) ; 



siml_3 = (finalO_3M«l) + (guessl_2M«l ) + (guessl_4M«l ) + (gues s2_3M«l ) ; 

sim2_3 = guessl_2M+guess3_2M+guessl_4M+guess3_4M; 
sim2_3 += guess2_3M«2; // multiply by four 

sim3_3 = (guess2_3M«l) + (guess3_2M«l ) + (guess3_4M«l ) + (gues s4_3M«l ) ; 

// sim4 3 = guess3 21vl+guess5 2M+guess3 4M+guess5 _4M; invalid 
sim4_3 = guess3_2M+guess4_2Miguess3_4M+guess4_4M; 
sim4_3 +- guess4_3M<<2; // multiply by four 



int errl_3 - (imagel_3M«3) - siml_3 

int err2_3 = (image2_3M«3) - sim2_3 

int err3_3 = (image3_3M«3) - sim3_3 

int err4_3 = (image4_3M«3) - sim4_3 

int siml_2 = 0; // column 2 
int sim2_2 = 0; 
int sim3_2 = 0 ; 
int sim4 2=0; 



/ / column 3 
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if( isOddPixel ) 
{ 

siml_2 = f inalO_lM+guess2_lM+finalO_3M-fguess2_3M; 
siml_2 += guessl_2M«2; // multiply by four 

sim2_2 = (guessl_2M«l) + (guess2_lM«l) + (guess2_3M«l) + (guess3_2M«l) ; 

sim3_2 = guess2_lM+guess4_lM+guess2_3M+guess4_3M; 
sim3_2 += guess3_2M«2; // multiply by four 

//sim4 2 - (guess3_ 2M«i) + {guess4 1M«1 ) + (guess4_ 3M«1 ) + (guessS 2M«1) ; invalid 

sim4_2 = (guess3_2M«l) + (guess4_lM«l ) + (guess4 3M«1 ) + (guess 4 2M«1) ; 

} " " 

else 

{ 

siml_2 = (f inalO_2M«l) + (guessl_lM«l ) + (guessl_3M«l ) + (guess2_2M«l) ; 

sim2_2 = guessl_lM+guess3_lM+guessl_3M+guess3_3M; 
sim2_2 += guess2_2M«2; // multiply by four 

sim3_2 = (guess2_2M«l) + (guess3_lM«l ) + (guess3_3M«l ) + (guess4_2M«l) ; 

// sim4 2 » guess3 lM+guess5 __lM+guess3 3M+guess5 3M; invalid 
sim4_2 = guess3_lM+guess4_lM+guess3_3M+guess4_3M; 
sim4_2 += guess4_2M«2; // multiply by four 

} 

int errl_2 = (imagel_2M«3) - siml_2; // column 2 
int err2_2 = (image2__2M«3) - sim2_2; 
int err3_2 = (image3_2M«3) - sim3_2; 
int err4_2 = (image4_2M«3) - sim4_2; 

// compute guess column 3 

if ( isOddPixel ) 

{ 

temp = (errl_2»2) + (errl__3»l) + (errl_4»2) ; // 2x 4x 2x 

temp += (err2_2»l) + err2_3 + (err2_4»l) ; // 4x 8x 4x 

temp += (err3_2»2) + (err3_3»l) + (err3_4»2) ; // 2x 4x 2x 

temp = temp * alphalM; 

temp = temp » 7; // numerator is assumed to be 4 

temp = temp + guess 2_3M; 

guess2_3M = max (0 , min (temp, 255) ) ; // clip value 

// err4___3 - don't update this value - we don't realy have enough information 
// I will need to check if an aproximation is better than nothing though 

} 

else 
{ 

//temp - (errO 2»2) + (errO 3»lj + (errO 4»2) ; // 2x 4x 2x 
temp = (finalO_3M-127)«3; 77 2x 4x 2x 

temp += (errl_2»l) + errl_3 + (errl_4»l) ; // 4x 8x 4x 

temp += (err2_2»2) + (err2_3»l) + (err2_4»2); // 2x 4x 2x 

temp - temp * alphalM; 

temp = temp » 7; // numerator is assumed to be 4 

temp = temp + guessl_3M; 

guessl_3M = max (0, min (temp, 255) ) ; // clip value 



temp = (err2_2»2) + (err2_3»l) + (err2_4»2) 
temp += (err3_2»l) + err3_3 + (err3_4»l) 

temp += (err4_2»2) + (err4_3»l) + (err4_4»2) 



// 2x 4x 2x 
// 4x 8x 4x 
// 2x 4x 2x 
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temp = temp * alphalM; 

temp = temp » 7; // numerator is assumed to be 4 

temp — temp + guess3_3M; 

guess3_3M = max { 0 , min (temp, 255 )) ; // clip value 

} 

// compute aim column 1 
int siml_l =0; // column 1 
int sim2_l =0/ 
int sim3_l = 0; 
int sim4_l = 0; 

if ( ! isOddPixel ) 
{ 



siml_ 
siml_ 


1 
1 


+= 


final0_0M+guess2_0M+final0_2M+guess2_2M; 
= guessl_lM«2; // multiply by four 




sim2_ 


1 




(guessl_lM«l) + (guess2_0M«l ) + (guess2_2M«l ) + (guess3_ 


_1M«1) ; 


sim3 
sim3_ 


1 
_1 


+= 


guess2_0M+guess4_0M+guess2__2M+guess4 2M; 
= guess3_lM«2; // multiply by four 




sim4 


1 




(guess3__lM«l) + (guess4_0M«l) + (guess4 2M«1) + (guess4 


_1M«1) ; 


siml_ 


_1 




(finalO_lM«l) + (guessl_0M«l) + (guessl_2M«l ) + (guess2_ 


_1M«1) ; 


sim2_ 
sim2 


1 
~1 


+= 


gue s s l_0M+gues s 3_0M+gues s l_2M+gues s 3 2M; 
= guess2_lM«2; // multiply by four 




sim3_ 


_1 




(guess2_lM«l) + (guess3_0M«l ) + (guess3_2M«l ) + (guess4 


_1M«1) ; 


sim4 
sim4_ 


1 
"l 


= gues s 3_0M+gue s s 4_0M+gue s s 3_2M+gue s s 4_2M; 
+= guess4_lM«2; // multiply by four 





} 



int errl_l = (imagel__lM«3) - siml_l; // column 1 

int err2_l = (image2_lM«3 ) - sim2_l; 

int err3_l = (image3_lM«3) - sim3_l; 

int err4_l = (image4_lM«3) - sim4_l; 

// Compute Guess for Column 2 
if( !isOddPixel ) 

{ "A 

temp = (errl_l»2) + (errl_2»l) + (errl_3»2) ; // 2x 4x 2x 

temp += (err2_l»l) + err2_2 + (err2_3»l) ; // 4x 8x 4x 

temp += (err3_l»2) + (err3 2»1) + (err3 3»2) ; // 2x 4x 2x 



temp = temp * alpha2M; 

temp = temp » 7; // numerator is assumed to be 4 

temp - temp + guess2_2M; 

gue s s 2_2M = max ( 0 , min ( temp ,255)); / / cl ip value 

// err4_3 ^ don't update this value - we don't realy have enough information 
// I will need to check if an aproximation is better than .nothing though 

} 

else 
{ 

// temp - (er.r0__l»2) + (err0_2»l) + (errO 3»2) ; // 2x 4x 2x 
temp = (final0_2M-127) «3; // 2x 4x 2x 

temp += (errl_l»l) + errl_2 + (errl_3»l) ; // 4x 8x 4x 



} 
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temp += (err2_l»2) + (err2_2»l) + (err2_3»2) ; // 2x 4jc 2x 
temp - temp * alpha2M; 

temp - temp » 7; // numerator is assumed to be 4 

temp - temp + guessl_2M; 

guessl_2M = max ( 0 , min (temp, 255 )) ; // clip value 

temp = (err2_l»2) + <err2_2»l) + (err2_3»2) ; // 2x 4x 2x 
temp += (err3_l»l) + err3_2 + (err3_3»l) ; // 4x 8x 4x 

temp += (err4_l»2) + (err4_2»l) + (err4_3»2) ; // 2x 4x 2x 

temp - temp * alpha2M; 

temp = temp » 7/ // numerator is assumed to be 4 

temp = temp + guess3_2M; 

guess3_2M = max (0, min (temp, 255) ) ; // clip value 



// Compute Guess for Column 1 

if ( isOddPixel ) 

{ 

temp = (errl_0M»2) + (errl_l»l) + (errl_2»2) ; // 2x 4x 2x 

temp += (err2_0M»l) + err2_l + (err2_2»l) ; // 4x 8x 4x 

temp += (err3_0M»2) + (err3_l»l) + (err3_2»2) ; /./ 2x 4x 2x 

temp = temp * alpha3M; 

temp = temp » 7; // numerator is assumed to be 4 

temp = temp + guess2__lM; 

guess2_lM = max (0, min (temp, 255) ) ; // clip value 

// err4 1 - don't update this^ value - we don't reaiy have enough information 



} 

else 
{ 



/'/ I will need to check if an aproximation is better than nothing though 



//temp - (errO 0M»2) + (errO + (errO 2»2) ; // 2x 4x 2x 

temp = (finalO_lM-127) «3; // 2x 4x 2x 

temp += (errl_0M»l) + errl_l + (errl__2»l) ; // 4x 8x 4x 

temp += (err2_0M»2) + (err2_l»l) + (err2_2»2) ; // 2x 4x 2x 

temp = temp * alpha3M; 

temp = temp » 7; // numerator is assumed to be 4 

temp = temp + guessl_lM; 

guessl_lM = max (0, min (temp, 255) ) ; // clip value 

temp = (err2_0M»2) + (err2_l»l) + (err2_2»2) ; // 2x 4x 2x 
temp += (err3_0M»l) + err3_l + (err3_2»l) ; // 4x 8x 4x 

temp 4-= (err4_0M»2) + (err4_l»l) + (err4_2»2); // 2x 4x 2x 

temp = temp * alpha 3M; 

temp = temp » 7; // numerator is assumed to be 4 

temp = temp + guess 3_1M; 

guess3_lM = max (0, min (temp, 255 )) ; // clip value 

} 

unsigned char rv = 0; 

if ( isOddPixel ) 
{ 

rv = guessl_0M; 

} 

else 
{ 

temp = (old_errl_0M»3) + (errl_0M»2) + (errl_l»3) ; //lx 2x Ix 
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temp = temp» 2; // divide by 4 

// make the signed, value fit in one byte I 

// I could do a. non-linear lookup table here! 

if ( temp < -127 ) temp = -127; 

if ( temp > 127 ) temp = 127; 

temp += 127; 

rv = (unsigned char) temp; 

} 

// Shift Error Values 
old_errl_0M = errl_0M; 
errl_0M = errl_l; 
err2_0M = err2_l; 
err3_0M = err3_l; 
err4_0M = err4_l; 

return rv; //return error value or subframe value 



